Exhaled Volatile Organic Compounds for Asthma Control Classification in Children with Moderate to Severe Asthma: Results from the SysPharmPediA Study.

RATIONALE: Early identification of children with poorly controlled asthma is imperative for optimizing treatment strategies. The analysis of exhaled volatile organic compounds (VOCs) is an emerging approach to identify prognostic and diagnostic biomarkers in pediatric asthma. OBJECTIVES: To assess the accuracy of gas chromatography-mass spectrometry based exhaled metabolite analysis to differentiate between controlled and uncontrolled pediatric asthma. METHODS: This study encompassed a discovery (SysPharmPediA) and validation phase (U-BIOPRED, PANDA). Firstly, exhaled VOCs that discriminated asthma control levels were identified. Subsequently, outcomes were validated in two independent cohorts. Patients were classified as controlled or uncontrolled, based on asthma control test scores and number of severe attacks in the past year. Additionally, potential of VOCs in predicting two or more future severe asthma attacks in SysPharmPediA was evaluated. MEASUREMENTS AND MAIN RESULTS: Complete data were available for 196 children (SysPharmPediA=100, U-BIOPRED=49, PANDA=47). In SysPharmPediA, after randomly splitting the population into training (n=51) and test sets (n=49), three compounds (acetophenone, ethylbenzene, and styrene) distinguished between uncontrolled and controlled asthmatics. The area under the receiver operating characteristic curve (AUROCC) for training and test sets were respectively: 0.83 (95% CI: 0.65-1.00) and 0.77 (95% CI: 0.58-0.96). Combinations of these VOCs resulted in AUROCCs of 0.74 ±0.06 (UBIOPRED) and 0.68 ±0.05 (PANDA). Attacks prediction tests, resulted in AUROCCs of 0.71 (95% CI 0.51-0.91) and 0.71 (95% CI 0.52-0.90) for training and test sets. CONCLUSIONS: Exhaled metabolites analysis might enable asthma control classification in children. This should stimulate further development of exhaled metabolites-based point-of-care tests in asthma.

Epigenomic response to albuterol treatment in asthma-relevant airway epithelial cells.

BACKGROUND: Albuterol is the first-line asthma medication used in diverse populations. Although DNA methylation (DNAm) is an epigenetic mechanism involved in asthma and bronchodilator drug response (BDR), no study has assessed whether albuterol could induce changes in the airway epithelial methylome. We aimed to characterize albuterol-induced DNAm changes in airway epithelial cells, and assess potential functional consequences and the influence of genetic variation and asthma-related clinical variables. RESULTS: We followed a discovery and validation study design to characterize albuterol-induced DNAm changes in paired airway epithelial cultures stimulated in vitro with albuterol. In the discovery phase, an epigenome-wide association study using paired nasal epithelial cultures from Puerto Rican children (n = 97) identified 22 CpGs genome-wide associated with repeated-use albuterol treatment (p < 9 × 10-8). Albuterol predominantly induced a hypomethylation effect on CpGs captured by the EPIC array across the genome (probability of hypomethylation: 76%, p value = 3.3 × 10-5). DNAm changes on the CpGs cg23032799 (CREB3L1), cg00483640 (MYLK4-LINC01600), and cg05673431 (KSR1) were validated in nasal epithelia from 10 independent donors (false discovery rate [FDR] < 0.05). The effect on the CpG cg23032799 (CREB3L1) was cross-tissue validated in bronchial epithelial cells at nominal level (p = 0.030). DNAm changes in these three CpGs were shown to be influenced by three independent genetic variants (FDR < 0.05). In silico analyses showed these polymorphisms regulated gene expression of nearby genes in lungs and/or fibroblasts including KSR1 and LINC01600 (6.30 × 10-14 ≤ p ≤ 6.60 × 10-5). Additionally, hypomethylation at the CpGs cg10290200 (FLNC) and cg05673431 (KSR1) was associated with increased gene expression of the genes where they are located (FDR < 0.05). Furthermore, while the epigenetic effect of albuterol was independent of the asthma status, severity, and use of medication, BDR was nominally associated with the effect on the CpG cg23032799 (CREB3L1) (p = 0.004). Gene-set enrichment analyses revealed that epigenomic modifications of albuterol could participate in asthma-relevant processes (e.g., IL-2, TNF-α, and NF-κB signaling pathways). Finally, nine differentially methylated regions were associated with albuterol treatment, including CREB3L1, MYLK4, and KSR1 (adjusted p value < 0.05). CONCLUSIONS: This study revealed evidence of epigenetic modifications induced by albuterol in the mucociliary airway epithelium. The epigenomic response induced by albuterol might have potential clinical implications by affecting biological pathways relevant to asthma.

Novel insights into the whole-blood DNA methylome of asthma in ethnically diverse children and youth.

BACKGROUND: The epigenetic mechanisms of asthma remain largely understudied in African Americans and Hispanics/Latinos, two populations disproportionately affected by asthma. We aimed to identify markers, regions and processes with differential patterns of DNA methylation (DNAm) in whole blood by asthma status in ethnically diverse children and youth, and to assess their functional consequences. METHODS: DNAm levels were profiled with the Infinium MethylationEPIC or HumanMethylation450 BeadChip arrays among 1226 African Americans or Hispanics/Latinos and assessed for differential methylation per asthma status at the CpG and region (differentially methylated region (DMR)) level. Novel associations were validated in blood and/or nasal epithelium from ethnically diverse children and youth. The functional and biological implications of the markers identified were investigated by combining epigenomics with transcriptomics from study participants. RESULTS: 128 CpGs and 196 DMRs were differentially methylated after multiple testing corrections, including 92.3% and 92.8% novel associations, respectively. 41 CpGs were replicated in other Hispanics/Latinos, prioritising cg17647904 (NCOR2) and cg16412914 (AXIN1) as asthma DNAm markers. Significant DNAm markers were enriched in previous associations for asthma, fractional exhaled nitric oxide, bacterial infections, immune regulation or eosinophilia. Functional annotation highlighted epigenetically regulated gene networks involved in corticosteroid response, host defence and immune regulation. Several implicated genes are targets for approved or experimental drugs, including TNNC1 and NDUFA12. Many differentially methylated loci previously associated with asthma were validated in our study. CONCLUSIONS: We report novel whole-blood DNAm markers for asthma underlying key processes of the disease pathophysiology and confirm the transferability of previous asthma DNAm associations to ethnically diverse populations.

The genomic history of the indigenous people of the Canary Islands.

The indigenous population of the Canary Islands, which colonized the archipelago around the 3rd century CE, provides both a window into the past of North Africa and a unique model to explore the effects of insularity. We generate genome-wide data from 40 individuals from the seven islands, dated between the 3rd-16rd centuries CE. Along with components already present in Moroccan Neolithic populations, the Canarian natives show signatures related to Bronze Age expansions in Eurasia and trans-Saharan migrations. The lack of gene flow between islands and constant or decreasing effective population sizes suggest that populations were isolated. While some island populations maintained relatively high genetic diversity, with the only detected bottleneck coinciding with the colonization time, other islands with fewer natural resources show the effects of insularity and isolation. Finally, consistent genetic differentiation between eastern and western islands points to a more complex colonization process than previously thought.

Human genetics influences microbiome composition involved in asthma exacerbations despite inhaled corticosteroid treatment.

BACKGROUND: The upper-airway microbiome is involved in asthma exacerbations despite inhaled corticosteroid (ICS) treatment. Although human genetics regulates microbiome composition, its influence on asthma-related airway bacteria remains unknown. OBJECTIVE: We sought to identify genes and biological pathways regulating airway-microbiome traits involved in asthma exacerbations and ICS response. METHODS: Saliva, nasal, and pharyngeal samples from 257 European patients with asthma were analyzed. The association of 6,296,951 genetic variants with exacerbation-related microbiome traits despite ICS treatment was tested through microbiome genome-wide association studies. Variants with 1 × 10-4 

Alpha-1 antitrypsin deficiency and Pi*S and Pi*Z SERPINA1 variants are associated with asthma exacerbations.

INTRODUCTION AND OBJECTIVES: Asthma is a chronic inflammatory disease of the airways. Asthma patients may experience potentially life-threatening episodic flare-ups, known as exacerbations, which may significantly contribute to the asthma burden. The Pi*S and Pi*Z variants of the SERPINA1 gene, which usually involve alpha-1 antitrypsin (AAT) deficiency, had previously been associated with asthma. The link between AAT deficiency and asthma might be represented by the elastase/antielastase imbalance. However, their role in asthma exacerbations remains unknown. Our objective was to assess whether SERPINA1 genetic variants and reduced AAT protein levels are associated with asthma exacerbations. MATERIALS AND METHODS: In the discovery analysis, SERPINA1 Pi*S and Pi*Z variants and serum AAT levels were analyzed in 369 subjects from La Palma (Canary Islands, Spain). As replication, genomic data from two studies focused on 525 Spaniards and publicly available data from UK Biobank, FinnGen, and GWAS Catalog (Open Targets Genetics) were analyzed. The associations between SERPINA1 Pi*S and Pi*Z variants and AAT deficiency with asthma exacerbations were analyzed with logistic regression models, including age, sex, and genotype principal components as covariates. RESULTS: In the discovery, a significant association with asthma exacerbations was found for both Pi*S (odds ratio [OR]=2.38, 95% confidence interval [CI]= 1.40-4.04, p-value=0.001) and Pi*Z (OR=3.49, 95%CI=1.55-7.85, p-value=0.003)Likewise, AAT deficiency was associated with a higher risk for asthma exacerbations (OR=5.18, 95%CI=1.58-16.92, p-value=0.007) as well as AAT protein levels (OR= 0.72, 95%CI=0.57-0.91, p-value=0.005). The Pi*Z association with exacerbations was replicated in samples from Spaniards with two generations of Canary Islander origin (OR=3.79, p-value=0.028), and a significant association with asthma hospitalizations was found in the Finnish population (OR=1.12, p-value=0.007). CONCLUSIONS: AAT deficiency could be a potential therapeutic target for asthma exacerbations in specific populations.

Epigenome-Wide Association Studies of the Fractional Exhaled Nitric Oxide and Bronchodilator Drug Response in Moderate-to-Severe Pediatric Asthma.

Asthma is the most prevalent pediatric chronic disease. Bronchodilator drug response (BDR) and fractional exhaled nitric oxide (FeNO) are clinical biomarkers of asthma. Although DNA methylation (DNAm) contributes to asthma pathogenesis, the influence of DNAm on BDR and FeNO is scarcely investigated. This study aims to identify DNAm markers in whole blood associated either with BDR or FeNO in pediatric asthma. We analyzed 121 samples from children with moderate-to-severe asthma. The association of genome-wide DNAm with BDR and FeNO has been assessed using regression models, adjusting for age, sex, ancestry, and tissue heterogeneity. Cross-tissue validation was assessed in 50 nasal samples. Differentially methylated regions (DMRs) and enrichment in traits and biological pathways were assessed. A false discovery rate (FDR) < 0.1 and a genome-wide significance threshold of p < 9 × 10-8 were used to control for false-positive results. The CpG cg12835256 (PLA2G12A) was genome-wide associated with FeNO in blood samples (coefficient= -0.015, p = 2.53 × 10-9) and nominally associated in nasal samples (coefficient = -0.015, p = 0.045). Additionally, three CpGs were suggestively associated with BDR (FDR < 0.1). We identified 12 and four DMRs associated with FeNO and BDR (FDR < 0.05), respectively. An enrichment in allergic and inflammatory processes, smoking, and aging was observed. We reported novel associations of DNAm markers associated with BDR and FeNO enriched in asthma-related processes.

Multi-omic approach associates blood methylome with bronchodilator drug response in pediatric asthma.

BACKGROUND: Albuterol is the drug most widely used as asthma treatment among African Americans despite having a lower bronchodilator drug response (BDR) than other populations. Although BDR is affected by gene and environmental factors, the influence of DNA methylation is unknown. OBJECTIVE: This study aimed to identify epigenetic markers in whole blood associated with BDR, study their functional consequences by multi-omic integration, and assess their clinical applicability in admixed populations with a high asthma burden. METHODS: We studied 414 children and young adults (8-21 years old) with asthma in a discovery and replication design. We performed an epigenome-wide association study on 221 African Americans and replicated the results on 193 Latinos. Functional consequences were assessed by integrating epigenomics with genomics, transcriptomics, and environmental exposure data. Machine learning was used to develop a panel of epigenetic markers to classify treatment response. RESULTS: We identified 5 differentially methylated regions and 2 CpGs genome-wide significantly associated with BDR in African Americans located in FGL2 (cg08241295, P = 6.8 × 10-9) and DNASE2 (cg15341340, P = 7.8 × 10-8), which were regulated by genetic variation and/or associated with gene expression of nearby genes (false discovery rate < 0.05). The CpG cg15341340 was replicated in Latinos (P = 3.5 × 10-3). Moreover, a panel of 70 CpGs showed good classification for those with response and nonresponse to albuterol therapy in African American and Latino children (area under the receiver operating characteristic curve for training, 0.99; for validation, 0.70-0.71). The DNA methylation model showed similar discrimination as clinical predictors (P > .05). CONCLUSIONS: We report novel associations of epigenetic markers with BDR in pediatric asthma and demonstrate for the first time the applicability of pharmacoepigenetics in precision medicine of respiratory diseases.

Classifying asthma control using salivary and fecal bacterial microbiome in children with moderate-to-severe asthma.

BACKGROUND: Uncontrolled asthma can lead to severe exacerbations and reduced quality of life. Research has shown that the microbiome may be linked with asthma characteristics; however, its association with asthma control has not been explored. We aimed to investigate whether the gastrointestinal microbiome can be used to discriminate between uncontrolled and controlled asthma in children. METHODS: 143 and 103 feces samples were obtained from 143 children with moderate-to-severe asthma aged 6 to 17 years from the SysPharmPediA study. Patients were classified as controlled or uncontrolled asthmatics, and their microbiome at species level was compared using global (alpha/beta) diversity, conventional differential abundance analysis (DAA, analysis of compositions of microbiomes with bias correction), and machine learning [Recursive Ensemble Feature Selection (REFS)]. RESULTS: Global diversity and DAA did not find significant differences between controlled and uncontrolled pediatric asthmatics. REFS detected a set of taxa, including Haemophilus and Veillonella, differentiating uncontrolled and controlled asthma with an average classification accuracy of 81% (saliva) and 86% (feces). These taxa showed enrichment in taxa previously associated with inflammatory diseases for both sampling compartments, and with COPD for the saliva samples. CONCLUSION: Controlled and uncontrolled children with asthma can be differentiated based on their gastrointestinal microbiome using machine learning, specifically REFS. Our results show an association between asthma control and the gastrointestinal microbiome. This suggests that the gastrointestinal microbiome may be a potential biomarker for treatment responsiveness and thereby help to improve asthma control in children.

Admixture mapping of severe asthma exacerbations in Hispanic/Latino children and youth.

BACKGROUND: In the USA, genetically admixed populations have the highest asthma prevalence and severe asthma exacerbations rates. This could be explained not only by environmental factors but also by genetic variants that exert ethnic-specific effects. However, no admixture mapping has been performed for severe asthma exacerbations. OBJECTIVE: We sought to identify genetic variants associated with severe asthma exacerbations in Hispanic/Latino subgroups by means of admixture mapping analyses and fine mapping, and to assess their transferability to other populations and potential functional roles. METHODS: We performed an admixture mapping in 1124 Puerto Rican and 625 Mexican American children with asthma. Fine-mapping of the significant peaks was performed via allelic testing of common and rare variants. We performed replication across Hispanic/Latino subgroups, and the transferability to non-Hispanic/Latino populations was assessed in 1001 African Americans, 1250 Singaporeans and 941 Europeans with asthma. The effects of the variants on gene expression and DNA methylation from whole blood were also evaluated in participants with asthma and in silico with data obtained through public databases. RESULTS: Genomewide significant associations of Indigenous American ancestry with severe asthma exacerbations were found at 5q32 in Mexican Americans as well as at 13q13-q13.2 and 3p13 in Puerto Ricans. The single nucleotide polymorphism (SNP) rs1144986 (C5orf46) showed consistent effects for severe asthma exacerbations across Hispanic/Latino subgroups, but it was not validated in non-Hispanics/Latinos. This SNP was associated with DPYSL3 DNA methylation and SCGB3A2 gene expression levels. CONCLUSIONS: Admixture mapping study of asthma exacerbations revealed a novel locus that exhibited Hispanic/Latino-specific effects and regulated DPYSL3 and SCGB3A2.

Medication use in uncontrolled pediatric asthma: Results from the SysPharmPediA study.

BACKGROUND: Uncontrolled pediatric asthma has a large impact on patients and their caregivers. More insight into determinants of uncontrolled asthma is needed. We aim to compare treatment regimens, inhaler techniques, medication adherence and other characteristics of children with controlled and uncontrolled asthma in the: Systems Pharmacology approach to uncontrolled Paediatric Asthma (SysPharmPediA) study. MATERIAL AND METHODS: 145 children with moderate to severe doctor-diagnosed asthma (91 uncontrolled and 54 controlled) aged 6-17 years were enrolled in this multicountry, (Germany, Slovenia, Spain, and the Netherlands) observational, case-control study. The definition of uncontrolled asthma was based on asthma symptoms and/or exacerbations in the past year. Patient-reported adherence and clinician-reported medication use were assessed, as well as lung function and inhalation technique. A logistic regression model was fitted to assess determinants of uncontrolled pediatric asthma. RESULTS: Children in higher asthma treatment steps had a higher risk of uncontrolled asthma (OR (95%CI): 3.30 (1.56-7.19)). The risk of uncontrolled asthma was associated with a larger change in FEV1% predicted post and pre-salbutamol (OR (95%CI): 1.08 (1.02-1.15)). Adherence and inhaler techniques were not associated with risk of uncontrolled asthma in this population. CONCLUSION: This study showed that children with uncontrolled moderate-to-severe asthma were treated in higher treatment steps compared to their controlled peers, but still showed a higher reversibility response to salbutamol. Self-reported adherence and inhaler technique scores did not differ between controlled and uncontrolled asthmatic children. Other determinants, such as environmental factors and differences in biological profiles, may influence the risk of uncontrolled asthma in this moderate to severe asthmatic population.

The upper-airway microbiome as a biomarker of asthma exacerbations despite inhaled corticosteroid treatment.

BACKGROUND: The response to inhaled corticosteroids (ICS) in asthma is affected by the interplay of several factors. Among these, the role of the upper-airway microbiome has been scarcely investigated. We aimed to evaluate the association between the salivary, pharyngeal, and nasal microbiome with asthma exacerbations despite receipt of ICS. METHODS: Samples from 250 asthma patients from the Genomics and Metagenomics of Asthma Severity (GEMAS) study treated with ICS were analyzed. Control/case subjects were defined by the absence/presence of asthma exacerbations in the past 6 months despite being treated with ICS. The bacterial microbiota was profiled by sequencing the V3-V4 region of the 16S rRNA gene. Differences between groups were assessed by PERMANOVA and regression models adjusted for potential confounders. A false discovery rate (FDR) of 5% was used to correct for multiple comparisons. Classification models of asthma exacerbations despite ICS treatment were built with machine learning approaches based on clinical, genetic, and microbiome data. RESULTS: In nasal and saliva samples, case subjects had lower bacterial diversity (Richness, Shannon, and Faith indices) than control subjects (.007 ≤ P ≤ .037). Asthma exacerbations accounted for 8% to 9% of the interindividual variation of the salivary and nasal microbiomes (.003 ≤ P ≤ .046). Three, 4, and 11 bacterial genera from the salivary, pharyngeal, and nasal microbiomes were differentially abundant between groups (4.09 × 10-12 ≤ FDR ≤ 0.047). Integrating clinical, genetic, and microbiome data showed good discrimination for the development of asthma exacerbations despite receipt of ICS (AUCtraining: 0.82 and AUCvalidation: 0.77). CONCLUSION: The diversity and composition of the upper-airway microbiome are associated with asthma exacerbations despite ICS treatment. The salivary microbiome has a potential application as a biomarker of asthma exacerbations despite receipt of ICS.

Epigenome-wide association study of lung function in Latino children and youth with asthma.

INTRODUCTION: DNA methylation studies have associated methylation levels at different CpG sites or genomic regions with lung function. Moreover, genetic ancestry has been associated with lung function in Latinos. However, no epigenome-wide association study (EWAS) of lung function has been performed in this population. Here, we aimed to identify DNA methylation patterns associated with lung function in pediatric asthma among Latinos. RESULTS: We conducted an EWAS in whole blood from 250 Puerto Rican and 148 Mexican American children and young adults with asthma. A total of five CpGs exceeded the genome-wide significance threshold of p = 1.17 × 10-7 in the combined analyses from Puerto Ricans and Mexican Americans: cg06035600 (MAP3K6, p = 6.13 × 10-8) showed significant association with pre-bronchodilator Tiffeneau-Pinelli index, the probes cg00914963 (TBC1D16, p = 1.04 × 10-7), cg16405908 (MRGPRE, p = 2.05 × 10-8), and cg07428101 (MUC2, p = 5.02 × 10-9) were associated with post-bronchodilator forced vital capacity (FVC), and cg20515679 (KCNJ6) with post-bronchodilator Tiffeneau-Pinelli index (p = 1.13 × 10-8). However, these markers did not show significant associations in publicly available data from Europeans (p > 0.05). A methylation quantitative trait loci analysis revealed that methylation levels at these CpG sites were regulated by genetic variation in Latinos and the Biobank-based Integrative Omics Studies (BIOS) consortium. Additionally, two differentially methylated regions in REXOC and AURKC were associated with pre-bronchodilator Tiffeneau-Pinelli index (adjusted p < 0.05) in Puerto Ricans and Mexican Americans. Moreover, we replicated some of the previous differentially methylated signals associated with lung function in non-Latino populations. CONCLUSIONS: We replicated previous associations of epigenetic markers with lung function in whole blood and identified novel population-specific associations shared among Latino subgroups.

Multiomics analysis identifies BIRC3 as a novel glucocorticoid response-associated gene.

BACKGROUND: Inhaled corticosteroid (ICS) response among patients with asthma is influenced by genetics, but biologically actionable insights based on associations have not been found. Various glucocorticoid response omics data sets are available to interrogate their biological effects. OBJECTIVE: We sought to identify functionally relevant ICS-response genetic associations by integrating complementary multiomics data sets. METHODS: Variants with P values less than 10-4 from a previous ICS-response genome-wide association study were reranked on the basis of integrative scores determined from (1) glucocorticoid receptor- and (2) RNA polymerase II-binding regions inferred from ChIP-Seq data for 3 airway cell types, (3) glucocorticoid response element motifs, (4) differentially expressed genes in response to glucocorticoid exposure according to 20 transcriptomic data sets, and (5) expression quantitative trait loci from GTEx. Candidate variants were tested for association with ICS response and asthma in 6 independent studies. RESULTS: Four variants had significant (q value < 0.05) multiomics integrative scores. These variants were in a locus consisting of 52 variants in high linkage disequilibrium (r2 ≥ 0.8) near glucocorticoid receptor-binding sites by the gene BIRC3. Variants were also BIRC3 expression quantitative trait loci in lung, and 2 were within/near putative glucocorticoid response element motifs. BIRC3 had increased RNA polymerase II occupancy and gene expression, with glucocorticoid exposure in 2 ChIP-Seq and 13 transcriptomic data sets. Some BIRC3 variants in the 52-variant locus were associated (P < .05) with ICS response in 3 independent studies and others with asthma in 1 study. CONCLUSIONS: BIRC3 should be prioritized for further functional studies of ICS response.

Pharmacogenetic studies of long-acting beta agonist and inhaled corticosteroid responsiveness in randomised controlled trials of individuals of African descent with asthma.

BACKGROUND: Pharmacogenetic studies in asthma cohorts, primarily made up of White people of European descent, have identified loci associated with response to inhaled beta agonists and corticosteroids (ICSs). Differences exist in how individuals from different ancestral backgrounds respond to long-acting beta agonist (LABA) and ICSs. Therefore, we sought to understand the pharmacogenetic mechanisms regulating therapeutic responsiveness in individuals of African descent. METHODS: We did ancestry-based pharmacogenetic studies of children (aged 5-11 years) and adolescents and adults (aged 12-69 years) from the Best African Response to Drug (BARD) trials, in which participants with asthma uncontrolled with low-dose ICS (fluticasone propionate 50 µg in children, 100 µg in adolescents and adults) received different step-up combination therapies. The hierarchal composite outcome of pairwise superior responsiveness in BARD was based on asthma exacerbations, a 31-day difference in annualised asthma-control days, or a 5% difference in percentage predicted FEV1. We did whole-genome admixture mapping of 15 159 ancestral segments within 312 independent regions, stratified by the two age groups. The two co-primary outcome comparisons were the step up from low-dose ICS to the quintuple dose of ICS (5 × ICS: 250 µg twice daily in children and 500 µg twice daily in adolescents and adults) versus double dose (2-2·5 × ICS: 100 µg twice daily in children, 250 µg twice daily in adolescents and adults), and 5 × ICS versus 100 µg fluticasone plus a LABA (salmeterol 50 µg twice daily). We used a genome-wide significance threshold of p<1·6 × 10-4, and tested for replication using independent cohorts of individuals of African descent with asthma. FINDINGS: We included 249 unrelated children and 267 unrelated adolescents and adults in the BARD pharmacogenetic analysis. In children, we identified a significant admixture mapping peak for superior responsiveness to 5 × ICS versus 100 µg fluticasone plus salmeterol on chromosome 12 (odds ratio [ORlocal African] 3·95, 95% CI 2·02-7·72, p=6·1 × 10-5) fine mapped to a locus adjacent to RNFT2 and NOS1 (rs73399224, ORallele dose 0·17, 95% CI 0·07-0·42, p=8·4 × 10-5). In adolescents and adults, we identified a peak for superior responsiveness to 5 × ICS versus 2·5 × ICS on chromosome 22 (ORlocal African 3·35, 1·98-5·67, p=6·8 × 10-6) containing a locus adjacent to TPST2 (rs5752429, ORallele dose 0·21, 0·09-0·52, p=5·7 × 10-4). We replicated rs5752429 and nominally replicated rs73399224 in independent African American cohorts. INTERPRETATION: BARD is the first genome-wide pharmacogenetic study of LABA and ICS response in clinical trials of individuals of African descent to detect and replicate genome-wide significant loci. Admixture mapping of the composite BARD trial outcome enabled the identification of novel pharmacogenetic variation accounting for differential therapeutic responses in people of African descent with asthma. FUNDING: National Institutes of Health, National Heart, Lung, and Blood Institute.

Role of Sex on the Genetic Susceptibility to Childhood Asthma in Latinos and African Americans.

Asthma is a respiratory disease whose prevalence changes throughout the lifespan and differs by sex, being more prevalent in males during childhood and in females after puberty. In this study, we assessed the influence of sex on the genetic susceptibility to childhood asthma in admixed populations. Sex-interaction and sex-stratified genome-wide association studies (GWAS) were performed in 4291 Latinos and 1730 African Americans separately, and results were meta-analyzed. Genome-wide (p ≤ 9.35 × 10-8) and suggestive (p ≤ 1.87 × 10-6) population-specific significance thresholds were calculated based on 1000 Genomes Project data. Additionally, protein quantitative trait locus (pQTL) information was gathered for the suggestively associated variants, and enrichment analyses of the proteins identified were carried out. Four independent loci showed interaction with sex at a suggestive level. The stratified GWAS highlighted the 17q12-21 asthma locus as a contributor to asthma susceptibility in both sexes but reached genome-wide significance only in females (p-females < 9.2 × 10-8; p-males < 1.25 × 10-2). Conversely, genetic variants upstream of ligand-dependent nuclear receptor corepressor-like gene (LCORL), previously involved in height determination and spermatogenesis, were associated with asthma only in males (minimum p = 5.31 × 10-8 for rs4593128). Enrichment analyses revealed an overrepresentation of processes related to the immune system and highlighted differences between sexes. In conclusion, we identified sex-specific polymorphisms that could contribute to the differences in the prevalence of childhood asthma between males and females.

LTA4H rs2660845 association with montelukast response in early and late-onset asthma.

Leukotrienes play a central pathophysiological role in both paediatric and adult asthma. However, 35% to 78% of asthmatics do not respond to leukotriene inhibitors. In this study we tested the role of the LTA4H regulatory variant rs2660845 and age of asthma onset in response to montelukast in ethnically diverse populations. We identified and genotyped 3,594 asthma patients treated with montelukast (2,514 late-onset and 1,080 early-onset) from seven cohorts (UKBiobank, GoSHARE, BREATHE, Tayside RCT, PAGES, GALA II and SAGE). Individuals under montelukast treatment experiencing at least one exacerbation in a 12-month period were compared against individuals with no exacerbation, using logistic regression for each cohort and meta-analysis. While no significant association was found with European late-onset subjects, a meta-analysis of 523 early-onset individuals from European ancestry demonstrated the odds of experiencing asthma exacerbations by carriers of at least one G allele, despite montelukast treatment, were increased (odds-ratio = 2.92, 95%confidence interval (CI): 1.04-8.18, I2 = 62%, p = 0.0412) compared to those in the AA group. When meta-analysing with other ethnic groups, no significant increased risk of asthma exacerbations was found (OR = 1.60, 95% CI: 0.61-4.19, I2 = 85%, p = 0.342). Our study demonstrates that genetic variation in LTA4H, together with timing of asthma onset, may contribute to variability in montelukast response. European individuals with early-onset (≤18y) carrying at least one copy of rs2660845 have increased odd of exacerbation under montelukast treatment, presumably due to the up-regulation of LTA4H activity. These findings support a precision medicine approach for the treatment of asthma with montelukast.